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ABSTRACT 

Context. The X-ray emission from Herbig Be/Ae stars has so far defied an unequivocal explanation. In later-type T Tauri stars, X-rays are 
thought to be produced by magnetically trapped coronal plasma, although accretion-shock induced X-rays have also been suggested. In earlier- 
type (OB) stars, shocks in unstable winds are thought to produce X-rays. 

Aims. We present the first high-resolution X-ray spectrum of a prototypical Herbig star (AB Aurigae), measure and interpret various spectral 
features, and compare our results with model predictions. 

Methods. We use X-ray spectroscopy data from the XMM-Newton Reflection Grating Spectrometers and the EPIC instruments. The spectra 
are interpreted using thermal, optically thin emission models with variable element abundances and a photoelectric absorption component. We 
interpret line flux ratios in He-like triplet of O vn as a function of electron density and the UV radiation field. We use the nearby co-eval classical 
T Tauri star SU Aur as a comparison. 

Results. AB Aurigae reveals a soft X-ray spectrum, most plasma being concentrated at 1-6 MK. The He-like triplet reveal no signatures of 
increased densities as reported for some accreting T Tau stars in the previous literature. There are also no clear indications for strong abundance 
anomalies in the emitting plasma. The light curve displays modulated variability, with a period of » 42 hr. 

Conclusions. It is unlikely that a nearby, undetected lower-mass companion is the source of the X-rays. Accretion shocks close to the star should 
be irradiated by the photosphere, leading to alteration in the He-like triplet fluxes of O vn, which we do not measure. Also, no indications for 
high densities are found, although the mass accretion rate is presently unknown. Emission from wind shocks is unlikely, given the weak radiation 
pressure. A possible explanation would be a solar-like magnetic corona. Magnetically confined winds provide a very promising alternative. The 
X-ray period is indeed close to periods previously measured in optical lines from the wind. 
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1. Introduction 

Herbig Ae/Be stars, first defined bv lHerbidlll960l) . are young 
intermediate-mass (« 2-10 M ) stars predominantly lo- 
cated near star-forming dark clouds. They show emission 
lines in their optical spectra, and their placement in the 
HR diagram prove s that they are pre-main sequence stars 
IStrom et allll972h . Herbig stars may therefore be considered 
to be intermediate-mass analogs of low-mass T Tauri stars 
(TTS), and in particular of their accreting variant, the classi- 
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cal T Tau stars (CTTS). The analogy between Herbig stars and 
CTTS extends to infrared excess emission indicative of heated 
circumstellar material (e.g., disks) and photometric variability. 
Herbig stars are important in the study of star formation be- 
cause their evolutionary scenario is intermediate between that 
of low-mass and high-mass stars. In the former class, accretion 
phase, dispersal of the circumstellar disk, and onset of hydro- 
gen burning occur sequentially in time whereas in the latter, the 
stars enter their main-sequence phase while still being embed- 
ded and accreting, making their study much more difficult. 

Like their main-sequence descendants, B and A-type stars, 
Herbig stars are generally supposed to be radiative in their 
interiors. The lack of a thick convection zone would make 
the operation of a solar-like a - co dynamo impossible, and 
hence no magnetic fields are expected at the surface of these 
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stars except possibly fossil magnetic fields trapped in the star 
since the initial cloud contraction phase. However, transient 
convection may be present in some of these stars during a 
short phase of deuterium burning in a shell, opening the pos- 
sibilit y that some dynamo-g enerated, non-potential fields de- 
velop dPalla&Stahlerl fl993): alternatively, a dynamo powered 
by rotational shear energy may generate some surface mag- 
netic fields in rapidly rotating, accreting Herbig stars as well 
jTout & PringleL fl995l) . Although difficult to detect, magnetic 
fields have recently indeed been measured on several Herbig 
stars, with (l ongitudin al) field str e ngths up to a few 100 G 
jDonati et all fl99l iHubrig et all Eool IWade e"Hil Eool). 
In sp ectropolarimetry studies of 50 Herbig stars, I Wade et alJ 
(2005) detected magnetic fields in 5 of them and discussed 
these stars to be progenitors of the magnetic Ap/Bp stars. 
IPraderie et all i ll 9861) have suggested magnetic activity in AB 
Aur based on their observations of periodically variable blue 
wings in the Mg n line of this star. Whatever the generation 
mechanism of magnetic fields, it may be important to include 
circumstellar disks in the picture as the magnetic fields may 
attach to the inner border of the disk so that they directly in- 
teract with accreting material. Moreover, the presence of mag- 
netic fields may also be important in the presence of stellar 
winds, because the wind could then be magnetically confined 
and t he plasma could be shock- heated to X-ray temperatures 
(e.g., Babel & Montmerle, 1997). We will address this point in 
Sect.lO 

X-ray emission is among the best tracers of magnetic fields 
in stars, alth ough it equal l y wel l dia gnoses wind or accre- 
tion sh ocks. iDamiani et alJ i 19941) and IZinnecker & Preibischl 
i 19941) were the first to systematically study X-ray emission 
from Herbig stars with the Einstein X-ray observatory and 
ROS AT, respectively. They reported surprisingly high detection 
rates of 11/31 and 11/21, respectively. These (and subsequent) 
studies have investigated X-ray emission in the context of other 
stellar parameters, with the following principal results: i) The 
X-ray luminosity, Lx, increases with the effective temperature 
(T e ff) and the stellar bolometric luminosity (L»), although the 
ratio L x /L t =s 10~ 6 - 10~ 5 is higher than in O stars (L x /L t * 
10~ 7 ) but lower than in T Tau stars (L x /L, « 10~ 4 - 10~ 3 ). 
ii) Lx is not correlated with the projected equatorial veloc - 
ity v sin i ( Dam iani etali 1 19941 IZinnecker & Preibischl 1 1994 . 
thus pointing either to a magnetic dynamo saturation effect or 
to X-rays not related to dynamo-driven magnetic fields, iii) 
Lx does correlate with indicators of disks, accretion, and out- 
flow such as infrared excess, mass accretion rate M, and wind 
velocity or wind momentum flux. The latter finding may in- 
deed support a non-magnetic origin of the X-rays, for example 
shocks in unstable winds analogous to more massi ve O stars 
(Damiani et all 1 1 994 IZin necker & Preibischl Il994l) . 

On the other hand, variability, flares, and extremely high 
electron temp eratures > 1 MK clearl y favor magnet i c pro- 
cesses jHamaguchi et all l2000l l2005t ISkinner et all [2004; 
iGiardino et all l2004 IStelzer et all 12006a) whether based on 
dynamo-generated or fossil magnetic fields near the star, or 
star-disk magnetic fields. Surface convection plays an impor- 
tant role in transferring energy into magnetic fields by stirring 
the magnetic footpoints, a process that leads to coronal heat- 



ing and mass ejections on the Sun and in magnetically active 
stars. If the convection zones are shallow as in late-A or early 
F-type main-sequence stars, then the magnetic dynamos ap- 
pear to operate rather inefficiently, leading to modest X-ray 
luminosities o f cor onal sources that reveal very soft spectra 
JPanzera et all [l999). If convection is absent - as in main- 
sequence A-type stars - then fossil magnetic fields are unlikely 
to build up non-potential configurations although magnetic ac- 
tivity is widespread am ong chemically peculiar Bp/Ap stars 
(Drake et al., Il987l fl994). Alternatively, however, winding-up 
magnetic fields connecting the star with the inner circumstel- 
lar disk may episodically release energy through reconnection, 
thus heating plasma and possibly ejecting plasmoids th at con- 
tribute to jets often seen in young stars JHavashi et all 1 19961: 
IMontmerle etailbOOfilHamaguchi et alll2005l) . 

Alternatively, X-rays have been suggested to be formed in 
accretion shocks in the CTTS TW Hya based on its exception- 
ally soft X-ray spectrum, indications for high electron densi- 
ties (%; 10 13 cm" 3 ) , and anomalous abundances of N and Ne. 
ISwartz et alJ {2005 ) proposed a similar scenario for the Herbig 
star HD 163296 based on its unusual soft spectrum (kT * 0.5 
keV). 

A caveat in gener al is that the majority o f Herbig stars are 
binaries or multiples jFeigelson et allE003l) . A close, uniden- 
tified T Tau companion that hides in the strong optical light 
could easily produce the observed X-rays because X-ray lu- 
minosities and electron temperatures of Herbig stars are of- 
ten quite similar to those of T Tau stars (see for example 
Stel zer et alJl2006al) . This emission is commonly interpreted 
as solar-type coronal magnetic activity. The binary hypothe- 
sis has become a relevant m odel for flaring X-rays in at least 
two Herbig stars (MWC 297, Hamagu chi et alJ20 00: field stars 
discussed as alternative X-ray sources a nd companion discov- 
ered bv lVink et alJl2005t and V892 Tau. IGiardino et aljl2 004: 
an 1.5-2M c ompanion was discovered by Sm ith et alJ l2005). 
IStelzer et alJ J2006al) studied a sample of 17 Herbig stars ob- 
served with Chandra and detected 13 of them in X-rays. Of 
these 13 stars, 7 have a known visual or spectroscopic compan- 
ion that could be the source of the observed X-rays. Only 35% 
of the detected sources cannot be explained by known com- 
panions. Furthermore, in the latter work, X-ray properties of 
Herbig stars are found to be very similar to X-rays properties 
of CTTS, leading to two possibilities: either the mechanism of 
the X-ray generation is similar for the two types of stars, or the 
X-rays are generated by (partly unknown) low-mass compan- 
ions. 

We report here the first high-resolution X-ray spectrum of a 
Herbig star, namely AB Aurigae, obtained with the Reflection 
Grating Spectrometer (RGS) on board XMM-Newton. Only the 
high resolution spectrum gives access to He-like line triplets, 
which yield information on densities or UV radiation fields 
in X-ray emitting re gions. This, i n turn , constrains the source 
location and extent dBehar et all |2004J). Further, high resolu- 
tion spectroscopy enables accurate abundance determination 
using individual spectral lines (including lines of N and O 
forming at low temperatures and located at long wavelengths), 
and, together with CCD spectra obtained with the European 
Photon Imaging Cameras (EPIC), accurate thermal modeling. 
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Our analysis shows that AB Aur is another very soft source 
with a moderate X-ray luminosity, but there are no signs of 
increased electron densities nor elevated abundances. We will 
use the nearby, co-eval, classical T Tau star SU Aur as an ideal 
comparison star to pinpoint fundamental differences in these 
X-ray sources with largely differing interiors. 

The structure of our paper is as follows. First, we describe 
our target in Sect. |2] We then present our data reduction pro- 
cedures in Sect. [3] Sect. |4] presents our results from the X- 
ray spectroscopy, and Sect.[5]specifically analyzes information 
from the He-like triplet of O vn. We discuss possible models in 
Sect. [6] and conclude in Sect. [71 

2. AB Aur 

Table ^ summarizes the basic properties and the principal X- 
ray parameters of AB Aur. For comparison, the properties of 
the Herbig star HD 16329 6, that has been reported to reveal 
a soft X-ray spectrum ( Swart z et all 12004 . and parameters of 
the CTTS SU Aur are also reported (see also Sect.g]& Sect.|6j. 

New, preliminary estimates of the surface temperature 
(Jeff) and stellar luminosity (L„) of AB Aur have been de- 
rived using the method developed by Fitzpatrick & Massa 
Jl999h . By utilizing the available UV (IUE: SWP + LWP) 
through optical (photometric: UBV) data, the energy distribu- 
tion can be modeled with surface fluxes represent ed by R. L. 
Kurucz's ATLAS 13 atmospheric models dKuruczl[l993l) . See 
iDeWarfetalJ i2003) for a more detailed description of this 
procedure and how it was implemented for SU Aur. A com- 
plete description of this analysis, as it pertains to AB Aur, is in 
preparation (DeWarf et al., in preparation). 

The rotation period of AB Aur is controversia l. Recent ob- 
serva tions suggest an inclination angle i ~ 21.5° dCorder et alJ 
2005 and references therein), that would suggest, using vsim 
and radius from Tabled a rotation period P = 12.9 hr. For 
an extreme value of i - 70° previously reported for the disk 
inclination (see references in ICorder et al.ll2005l) . we obtain 
P — 33 hr. A period o f this order is supported by modulations in 
Mgn and Can lines JPraderie et al 1 fl986l ICatala et all ll986T) 
and in photospheric lines llCatala et all 1 19991) that reveal peri- 
ods P = 32 - 34 hr and P = 43 - 45 hr (see Sect. O for a 
detailed discussion). 

To our knowledge, no direct detection of surface magnetic 
fields of AB Aur has so far been obtained. ICatala et alJ (11999) 
estimated an upper limit to the strength of the magnetic field of 
300 G. 

3. Observations and data analysis 

AB Aur was detected in an XMM-Newton Jjansen et alll200ll) 
observation pointing at the nearby CTTS SU Aur (separation 
between SU Aur and AB Aur: « 2.5')- The observation was re- 
trieved from the archive as part of the XM M-Newton Extended 
Survey of Taurus (XEST) described in Giid el et all (2006) 
(XEST observation number of AB Aur: XEST-26-043, and of 
SU Aur: XEST-26-067). Table[2]summarizes observing param- 
eters. The two European Photon I maging Cameras (EPICs) of 
the MOS type jTurner et alll200ll) and the Reflection Grating 



Table 1. Parameters for AB Aur, SU Aur, and HD 163296 



Parameter 


AB Aur 


HD 163296 


SU Aur 


Magnetic field 


Y?(l) 


Y?(17) 


Y 


Spectrum 


B9.5e-A0 (2,3) 


AlVe (4) 


G2 (5) 


U [Lq] 


49 (16) 


30(4) 


9.9 (6) 


Mass [M G ] 


2.7(16) 


2.3 (4) 


1.9 (X) 


Radius [R a ] 


2.3 (16) 


2.1 (4) 


3.1 (X) 


A v [mag] 


0.25 (23) 


0.25 (4) 


0.5 (7) 




0.47 (16) 






Teff [K] 


9750 (8) 


9300 (4) 


5860 (6) 




10050 (16) 






Age [ Myr] 


4(16) 


4(9) 


4(7) 


Disk 


Y(21) 


450 AU (9) 


Y 


Companions? 


evidence (20) 


not det. (9) 




Radio spectrum 


wind (11) 


wind (12) 




EW(Hor) [A] 


27 (13) 




2-6 (X) 


vsini [km s _1 ] 


80 (18) 


120 (19) 


65 (22) 


kT [keV] 


0.46 (25) 


0.49 (9) 


1.9 (25) 


logL x [erg s -1 ] 


<30.3 (3) 


29.6 (9) 


30.9 (25) 




29.5 (10) 








29.6 (25) 






logLx/L„ 


-5.8(10) 


-5.48 (9) 


-3.7 (25) 




-5.6 (25) 






\ogN H [cm" 2 ] 


20.7 (24, X) 


20.88 (9) 


21.4 (25) 


X-ray variable? 


slow (25) 


marginal? (9) 


flares (24) 



REFERENCES given in parentheses: 1 Prad erie et alJ C 



- see text for further discussion; magnetic field assumed 


for CTTS. 2 Hamann & 


PerssorJ Jl992l): 3 iHamaauchi et alJ 


12005) (and references therein); 4 van den Ancker et al. 1 1998); 5 


Kenvon & Hartmann 1 1995; 


i: 6 iLuhmarJ <2004: 7 IDeWarfetalJ 


<2003l): 8lAcke & Waelkensf 


<2004l): 9lSwartz et alJ fe005l) (and ref- 



(1993); 12 Blondeletal. ( 199 31) 13 iHerb ig & Bell 11983) : 14 
iGradv et alJ fcOOOl) : 15 iDevine et alJ l200d): 16 deWar f & 



Fitzpatrick, private communication; 



17 iDeleuil et alJ (2005); 



18 



iBohm & Catai3jl993h: 19lHillenbrand et alJ<1992l): 20lBaines et alJ 
<2006l): 21 ICorder et al] |2005l): 22lHartmann & Stauffed <1998l) : 23 
Rober ge et alJfeOOll) : 24 Robrade & Schmit tH2006j): 25 t his work. X 
see XEST survey results and tableThl lGudeleHI 11 120061) for deriva- 
tions and references. 



Spectrometers (RGSs: lden Herder et al . 2001) were active dur- 
ing the observation, while the EPIC PN camera was out of oper- 
ation. Both MOS instruments observed in full frame mode and 
used the thick filter to suppress excessive optical load from AB 
Aur. The EPIC detectors operate in the energy range of 0.15- 
15.0 keV with a medium spectral resolution of approximately 
E/AE = 20-50. The RGSs are suited for high-resolution spec- 
troscopy in the wavelength range of 6-35 A and have a reso- 
lution of A/1 =s 60 - 76 mA. The RGS detectors contained the 
dispersed spectra of both AB Aur and SU Aur with sufficient 
separation to make mutual spectral contamination negligible 
(see belo w). The spectrum of SU A ur has recently been dis- 
cussed by Rob rade & Schmittl J2006I) . 

The data were reduced using the Science Analysis System 
(SAS) version 6.1. The EPIC MOS data were reduced us- 
ing the task emchain and the sources were detected using 
the maximum likelihood detection algorithm emldetect (see 
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Table 2. Observing log 



Table 3. Wavelength ranges used for the spectral fitting 



XMM-Newton ObsID 
Boresight RA (J2000.0) 
Boresight 6 (J2000.0) 
Start time (UT) 
Stop time (UT) 
Exposure time (s) 
MOS Mode and Filter 



0101440801 
4 h 55 h 59.0 s 
30° 34' 02" 
2001-09-21 01:34:17 
2001-09-22 13:34:31 
129614 

Full Frame, thick filter 



Idudel et aljEoOfil for further details). To extract the two spec- 
tra from each RGS detector, we proceeded as follows. We 
first applied the standard processing performed by the RGS 
metatask rgsproc to each source position. We extracted the to- 
tal (source+background) spectra and the background spectra 
separately. For both SU Aur and AB Aur, we included 85% of 
th e cross-dispersion Poin t Spread Function (PSF) as also done 
bv lTelleschi et alJ J2006h (xpsfincl = 85). For the background, 
the exclusion region with regard to the cross-dispersion PSF 
and the inclusion region of the pulse-height distribution were 
kept at default values (95% of the PSF and 90% of the pulse- 
height distribution, respectively, i.e., xpsfexcl = 95 and pdistincl 
= 90). We verified that this way the spatial extraction regions 
on the detector were adjacent to each other but just not over- 
lapping. This means that each spectrum collects approximately 
(100-85)/2 = 7.5% of the counts of the other spectrum. The 
contamination is such that counts from SU Aur are shifted in 
wavelength by approximately -0.4 A in the AB Aur spectrum 
within the wavelength range of interest here (12-22 A), and 
contaminating counts from AB Aur are shifted by +0.4 A in 
the SU Aur spectrum. 

The background defined for each source outside its source 
extraction region is now still contaminated by the other source. 
Therefore, we defined the coordinate of the secondary source 
(AB Aur for the SU Aur spectrum and vice versa), added them 
to the source list using the rgssources task, and computed a 
new extraction map, with the secondary source excluded from 
the background, using the rgsregion task. Finally, we extracted 
each spectrum again, using only the regions outside the SU Aur 
and AB Aur source regions for the background spectra. The 
background for the AB Aur spectrum is extracted outside the 
two adjacent source regions. The background portion on the far 
side of the SU Aur spectrum again comprises approximately 
7.5% of the SU Aur source counts. Because this background 
spectrum is subtracted from the total spectrum extracted at the 
AB Aur position, the contamination is approximately corrected 
for. Further considering the low S/N, the similar line fluxes in 
both spectra, and the line-dominated spectrum of AB Aur (see 
below, while most of the counts in SU Aur are in continuum 
and therefore distributed in wavelength), the mutual contami- 
nation is negligible and much below the noise level. For exam- 
ple, no significant mutual contamination is seen by the strong 
O vm lines above the noise level in either of the spectra (see 
presentation of the RGS spectra in Fig.[5]below). The periods 
affected by high background flaring were excluded from the 
spectral analysis. 



Instrument 


AB Aur 


SU Aur 


MOS1 




1.5-9.35 A 


MOS2 


1.5-9.35 A 




RGS1 


10.0-28.0 A 


8.3-25.0 A 


RGS2 


8.3-26.5 A 


8.3-25.0 A 



Because we wanted to put weight on the high-resolution 
spectra obtained with RGS, we performed the data analysis on 
both RGS spectra, adding only the short- wavelength portion of 
one of the two MOS detectors in order to access line features 
of Mg, S i and Fe that are im portan t for our abundance an aly- 
sis (see lTelleschi et al.lE o05). As in lTelleschi et alJ ll2006i) . we 
used MOS 1 for SU Aur for this purpose, while for AB Aur, we 
preferred MOS2 because an unidentified background feature 
distorted the MOS1 spectrum around a wavelength of 9.3 A 
(not present in MOS2). The wavelength intervals of each in- 
strument used for the data analysis are summarized in Table [5] 
As a check, we confirmed our results by using the entire useful 
energy range of MOS (0.2-10 keV), and the results are consis- 
tent with those reported here. 

We fitted the spectra in XSPEC jArnaudt 1 19961) us- 
ing optically-thin collisionally-ionized plasma models calcu- 
lated with the As trophysical Plasma Emission Code (APEC; 
ISmith et alJl200lh . In order to account for calibration discrep- 
ancies between the RGS and the MOS instruments, we added 
constants as effective-area renormalization factors. These fac- 
tors were fixed at 1.0 for MOS and at 1.05 for both RGS (see 
iKirsch et all2004 . 

Since we used the ^-statistic for our spectral fitting, we 
binned the spectra to a minimum of 20 counts per bin for the 
RGS and a minimum of 15 counts per bin for the MOS. The 
resulting bin width varies between 0.04 A (in the Ovn line) 
and 0.73 A (just longward of the Ovm Lyo- line) in the RGS, 
and between 0. 15 A and 1 .4 A in the MOS spectrum. 

For AB Aur the hydrogen column density Nh was fixed 
at the value found in the XEST survey analysis, namely 5 x 
10 20 cirT 2 , which is consistent with the recent me asurements 
of Ay — 0.25 mag reported by Roberg e et alJ(l200ll) . assuming 
a standard conversi on Nh ~ 2.0 x 10 21 Ay applicable to the 
interstellar medium JVuong et al.l2003l and references therein). 
This value agrees q uite well with the Nh found explicitly by 
iRoberge et all feOOll) . Nh = 4.4 x 10 20 cirT 2 . For SU Aur, we 
needed to fit Nh in order to get a good fit to the RGS spectra. 

The spectra were fitted with two different models: a model 
describing a differential emission measure distribution (EMD) 
approxim ated by two power laws as used in the XEST survey 
analysis (Gii del et all I2006I) . and a model with two or three 
isothermal plasma components. For both models we applied a 
procedure in which we simultaneously fitted large wavelength 
intervals of the three spectra with template spectra computed in 
XSPEC. Alternative, iterative methods based on extracted line 
fluxes are not feasible given the low S/N ratio of our spectra. 
There are only few explicitly measurable lines, and each set 
of lines from a given element is confined to narrow formation 
temperature interval. Also, such methods cannot be applied to 
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MOS CCD data. Previous studies have shown excellent agree- 
men t between iterativ e methods and the method applied here 
( lTelleschietalll2005h . The EMD model is approximated by a 
grid of isothermal components spaced regularly by 0. 1 dex in 
log T in such a way that the lower- T and the higher- T portions 
are each described by a power law. This model has been sug- 
gested from our previous work on high-resolutio n X-ra y spec- 
troscopy of young solar analogs (Telle schi et all I2Q05H . and is 
described by 



Q(T) = 



EM ■ (T/T r 
EM ■ (T/Taf 



JorT < T 
JorT> T 



0.035 
0.030 
0.025 
0.020 
0.015 
0.010 
0.005 

8:828 



0.015 - 



soft (0.3-1 keV) 



(1) S 

2 0.010 



hard (1-4 keV) 



where To is the temperature where the two power-laws cross, 
and EMo is the EM per log T at this crossing point. The slopes 
of the power laws below and above To are a and ft, respec- 
tively. The power laws have high-energy and low-energy cut- 
offs at log T = 8 and log T — 6, respectively. Given the 
poor energy resolution of CCDs below 1 keV, we had fixed 
a at a value of +2 throughout the XEST survey, and we will 
do so here for comparison reasons. A slope of 2 is compat- 
ible with findings from DEM analyses in other pre-main se- 
quence and main-seq uence active stars llArgiroffi et all [2004; 
iTelleschi et alll2005l) . 

However, because the X-ray emission of a Herbig Ae/Be 
star could be non-coronal, or the coronal thermal structure 
could be unlike that in late-type T Tau stars, it is possible that 
the simplified EMD model proposed here is inappropriate. In 
a separate approach, we will therefore also fit a. Further, be- 
cause the power-law EMD structure might be inappropriate 
altogether, we will test our results using a multi-temperature 
model that can reveal the temperatures where most EM re- 
sides. We therefore analyzed the AB Aur spectra with a 2- 
component model. A third component was not needed, given 
the low signal-to-noise ratio of the spectrum and, as we will de- 
scribe below, its rather narrow range of temperatures in which 
emitting plasma is found. Because spectral lines are formed 
over an interval of typically 0.3 dex in temperature, a temper- 
ature range of 2-7 MK (the range of formation temperatures 
present in the spectrum) will require no more than two ther- 
mal components. We also analyzed the SU Aur spectra with a 
3-component model. 

For both approaches, we fitted the abundances of the lines 
observed in the spectra simultaneously with the thermal mod- 
els. Abundances that do not show significant features in the 
spectra, such as those of C, S, Ar, Ca and Al (and N for 
SU Aur), were fixed at the values used in the general XEST 
survey data analysis (C =0.45, N=0 788, S= 0.417, Ar=0.55, 
Ca=0.195, and Al=0.5. ICiiidel et alJl2006l) . These abundances 
were thus arranged in such a way that they describe a weak 
"inverse First Ionization Potential Effect" (higher FIP implies 
higher abundances relative to the photospheric values, the latter 
assumed to be solar), as often observed in young stars 

Finally, we computed the X-ray luminosities, Lx, from the 
spectral model in t he energy range of 0. 3-10.0 keV, assuming a 
distance of 140 pc JPeWarf et alJ2003l and references therein). 
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Fig. 1. Light curves of AB Aur, combining counts from both 
MOS detectors, binned to 5 ks. Panels from top to bottom: 
Soft counts in the 0.3-1 keV range; hard counts in the 1-4 keV 
range; counts in the combined range, 0.3-4 keV; hardness, de- 
fined by the ratio of hard to soft counts. 



4. Results 

4.1. Light curves 

We first present the X-ray light curve of AB Aur in Fig. ^ 
We have co-added counts from both EPIC MOS detectors af- 
ter background subtraction. The background contribution to the 
total light curve was in fact negligible except during the final 
5 ks when a background flare occurred that reached » 50 % 
of the total count rate. The observed variability can therefore 
not be attributed to imperfect background treatment but is in- 
trinsic to the stellar source. We see slow variability on time 
scales of about one day, by somewhat less than a factor of two. 
In the variability an alysis of the XEST sources presented by 
IStelzer etaD d2006bl) . the variablity of AB Aur was confirmed 
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Fig. 2. EPIC spectra of SU Aur (MOS1) and AB Aur (MOS2). 
Important lines are labeled. 

by the Kolmogorov-Smirnov test. No significant trend is seen 
in the hardness, defined as the ratio between the count rates in 
the hard (1-4 keV) and the soft (0.3-1 keV) band. We have fit- 
ted the light curve with a sine function which is also plotted in 
Fig-E The fit is excellent (x 2 = 12.9 for 22 dof) with a period 
of 42.2^3 ^ hr in the total band (in the soft and hard band we 
find 48.5^'g 6 hr and 40.4^2 hr, respectively). A similar mod- 
ulation of 42 - 45 hr was found previously for the He i and 
Mg ii lines (see discussion below). 

We have studied the light curve variability with a statistical 
test, to assess the presence of a flare in the last 30 ks, where 
an increase in the count rates can be seen, particularly in the 
hard and total light curves (second and third panel in Fig. ^ 
respectively). We first tested our data against a constant count- 
rate model using the x 1 -statistic. We find the probability for 
this part of the light curve to be constant, P(const)=0.73 for 
the total band, and P(const)=0.30 for the hard band. As a sec- 
ond test, we used the Kolmogorov-Smirnov statistic, to obtain 
P(const)=0.09 and P(const)=0.07 for the total and hard light 
curves, respectively. The variability of the last 30 ks in the light 
curves is thus at best marginal. 

4.2. Spectra 

Fig. |2 compares the background-subtracted EPIC MOS spec- 
tra of AB Aur and SU Aur. A number of spectral differences 
are obvious. The spectrum of SU Aur reveals signatures of 
a very hot coronal plasma, with outstanding lines of Mgxi, 
Mg xn, Si xm and, most notably, the Fe K complex mostly due 
to Fe xxv at 6.7 keV. A steep drop toward the lowest energies 
indicates considerable photoelectric absorption. In contrast, the 
AB Aur spectrum falls off rapidly above 1 keV, showing its 
peak flux around 0.7-0.9 keV. This flux peak is mainly due to 
lines of Fe xvn. These spectral properties let us anticipate a soft 
source for AB Aur, and the shallow fall-off toward low energies 
indicates rather low photoelectric absorption. 

Fig-EJshows fluxed, co-added RGS1+2 spectra of SU Aur 
(top) and AB Aur (bottom) in the line dominated region. The 



spectra have been rebinned to a bin width of 0.042 A for AB 
Aur and 0.035 A for SU Aur. These spectra further corroborate 
the differences between the two X-ray sources. Whereas SU 
Aur reveals a strong continuum, indicating a hot source, there 
is little evidence for continuum emission in AB Aur. Further, 
the flux ratio between the Nex Lya line at 12.1 A (formation 
temperature 6.3 MK) and the Neix resonance line at 13.44 A 
(formation temperature 4 MK) is considerably higher in SU 
Aur than in AB Aur, again emphasizing the dominance of hot 
plasma in the former. Furthermore, the Neix line feature at 
13.5 A in SU Aur is dominated by Fexix, which is formed 
at higher temperatures (formation temperature 7.9 MK). Two 
further features are striking: First, the Fe lines of SU Aur are 
very strong, dominating the spectrum and comparing in flux 
with the O vm Lya and the Ne x Lya lines (although the O vm 
line is partly suppressed by photoelectric absorption). Such 
line ratios are unusual among very active, main-sequence solar 
analogs where Fe line fluxes are modest due to low Fe abun- 
dance dTelleschi et all E"005h . We thus anticipate an unusually 
high abundance of Fe in this spectrum. 

The second feature of interest here is the unusually high 
flux of the O vii lines of AB Aur, despite some photoelectric ab- 
sorption that suppresses the flux at these wavelengths. The total 
O vii flux appears to be similar to the flux in the O vm Lye line, 
which is a property of very inactive stellar coronae with tem- 
peratures of 3-5 MK ( Telles chi etaD.120051) . (Note that the fea- 
ture at 18.5 A is neither coincident with the O vii He/? line nor 
due to contamination by the O vm Lye line of SU Aur. Those 
two line features would both be slightly but significantly long- 
ward of this wavelength, namely at 18.6 A. The 18.5 A feature 
is due to a 3<x spike exclusively in the RGS2 detector and is 
therefore spurious). We note that the O vii lines are not present 
in the spectrum of SU Aur due to considerable photoelectric 
absorption. 

4.3. Thermal structure 

We now present the numerical results of our spectral fits to 
these data. Table |4] lists the numerical results for both mod- 
els and both stars. In Figure |4] we plot the data together with 
the best fit of the EMD model for AB Aur. 1 We define the av- 
erage coronal temperature, T- dv , as the logarithmic average of 
all temperatures used in the model, applying the correspond- 
ing EM as weights. This measure corresponds to the electron 
temperature itself for an isothermal plasma. For a continuous 
EMD as discussed here, r av represents a temperature grossly 
characteristic of the spectral shape. This is also true for multi- 
component plasma (e.g., 2-T or 3-T plasmas) although there 
might be no plasma present at T av itself. This also occurs when 
1-T fits are made to low-quality coronal spectra. 



1 The low bin at 16 A in RGS2, coincident with the O vm+Fe xvm 
lines, corresponds to a CCD gap. At the same time, somewhat in- 
creased flux in the RGS1 spectrum slightly longward of the two 
lines, probably due to imperfect background subtraction, increases the 
RGS1 flux in the respective bin in Fig.|4|above our fit. We checked at 
higher resolution that the O vm+Fe xvm line is indeed correctly fitted. 
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Fig. 3. Fluxed, coadded RGS1+RGS2 spectra of SU Aur (top) and AB Aur (bottom). The spectra are background subtracted. 
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Both the multi-temperature and the EMD models provide 
similar results: T in and the X-ray luminosities are similar. The 
2T or 3r-model fits and the EMD model fits show similar^ 2 
values, and are therefore of statistically equal quality. Further, 
most of the abundance values agree favorably between the two 
methods, confirming the robustness of our results. 



For AB Aur we report the results for EMD models with 
a frozen at 2 and with a free (second and third columns of 
Tabled respectively). The peak temperature Tq increases from 
4.4 MK for a fixed at 2 to 6.3 MK for free a. However, the main 
results are in very close agreement for the two approaches: the 
abundances are consistent within the error bars, the X-ray lu- 
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Table 4. Results of spectral model interpretation 1 



AB Aur SU Aur 





2T model 


DEM model 2 


DEM model 3 


3T model 


DEM model 2 


N H [10 21 cm- 2 ] 


= 0.5 2 


= 0.5 4 


= 0.5 4 


3.1 (3.0,3.3) 


3.2 (3.0, 3.4) 


N 5 


0.45 (0.22, 0.78) 


0.57 (0.30, 1.13) 


0.58 (0.27, 1.13) 


- 


- 


O 5 


0.20 (0.16, 0.42) 


0.22 (0.13,0.32) 


0.27 (0.17, 0.45) 


0.50 (0.40, 0.61) 


0.30 (0.24, 0.38) 


Ne 5 


0.60 (0.41, 0.88) 


0.62 (0.48, 1.04) 


0.77 (0.49, 1.22) 


0.61 (0.46, 0.78) 


0.38 (0.28, 0.53) 


Mg 5 


0.21 (0.09, 0.38) 


0.28 (0.13,0.74) 


0.34 (0.14, 0.62) 


1.51 (1.33, 1.74) 


1.17 (1.05, 1.43) 


Si 5 


0.70 (0.48, 1.01) 


0.90 (0.60, 1.32) 


0.96 (0.64, 1.51) 


0.79 (0.69, 0.92) 


0.64 (0.57, 0.79) 


S 5 


- 


- 


- 


0.64 (0.49, 0.79) 


0.57 (0.45, 0.72) 


Fe 5 


0.23 (0.17, 0.32) 


0.29 (0.22, 0.47) 


0.37 (0.24, 0.61) 


0.81 (0.73,0.91) 


0.67 (0.61, 0.77) 


Ti [MK] 


2.45 (2.10, 2.81) 


- 


- 


7.54 (7.30, 7.75) 


- 


T 2 [MK] 


6.99 (6.62, 7.41) 


- 


- 


13.97 (12.05, 18.06) 


- 


J 3 [MK] 


- 


- 


- 


38.14 (35.68,41.14) 


- 


EM, [10 52 cm" 3 ] 


2.11 (0.87, 3.67) 


- 


- 


11.27 (9.18, 14.04) 


- 


EM 2 [10 52 cm" 3 ] 


3.44 (3.06, 4.31) 


- 


— 


7.25 (4.53, 14.35) 


- 


EM 3 [10 52 cm" 3 ] 


- 


- 


- 


30.49 (24.66, 33.21) 


- 


To [MK] 


- 


4.38 (2.69, 5.68) 


6.33 (5.26, 8.32) 


- 


7.66 (6.95, 8.13) 


a 


- 


= 2.00 


0.47 (0.00, 1.20) 


- 


= 2.00 


P 




-1.9 (-2.57, -1.52) 


-2.05 (-3.00, -1.57) 




-0.05 (-0.11, 0.06) 


EM [10 52 cm- 3 ] 6 




5.21 


5.08 




57.84 


T av [MK] 


4.69 


4.71 


4.07 


22.65 


20.07 


L x [10 30 erg/s] 7 


0.40 


0.39 


0.39 


7.40 


7.79 


L x [10 30 erg/s] 8 


0.55 


0.53 


0.56 


8.53 


9.15 




1.06 


1.02 


1.00 


1.20 


1.23 


dof 


78 


79 


78 


443 


446 



68% error ranges are given in parentheses. 

2 Fitted with a held fixed at 2 like in the XEST survey (Giidel et al.. 2006') . 

3 Fitted with a being a free parameter. 

4 Held fixed at value found in the XEST survey (Giidel et al., 2006). 

5 Element abundances are with respect to solar values given by Anders & Grevesse (1989) (Grevesse & Sauval 1999 for Fe). 

6 Total EM integrated over temperature bins between log T = 6 - 7.9 [K] (see Giidel et al. 2006 for more details). 

7 Determined in the 0.3-10.0 keV band. 

8 Determined in the 0.1-10.0 keV band. 



minosities are essentially unchanged, and T av is only slightly 
smaller when we fit a. This suggests that the two solutions are 
equivalent. On the other hand, the errors are larger if a is a free 
fit parameter (except for the Mg abundance). Therefore, and 
for consistency with the fits of the EPIC spectra in the XEST 
survey, we henceforth use the results from the EMD fit with a 
fixed at a value of 2. 

The Nh value for SU Aur was found to be Nh = (3.1 - 
3.2) x 10 21 cm" 2 , w hich agrees well with the result reported by 
ISkinner& Walter I dl998l) from an ASCA observation, Nh - 
2.8 x 10 21 cm 2 an d is only somewhat h igher than expected 
from A v = 0.9 mag jKenvon & Hartmannlll995i) . 

The results confirm the peculiar thermal structure of the AB 
Aur source. The average temperature from the 2-T fit amounts 
to 4.7 MK, while the EMD model peaks at 4.4 MK and rapidly 
falls off toward higher temperatures, with a power-law slope of 
-1.9. Such low temperatures are unusual for coronae of young 
stars where usually temperatures in excess of 10 MK are found. 
An example is SU Aur: The 3-T model shows the largest 
amount of EM at the highest temperature at as 40 MK, with 
T av ss 23 MK. Again, the EMD model supports this finding, 
where we find Tq at 7.7 MK beyond which the EMD is nearly 



flat, resulting in an average temperature of 20 MK. Franciosini 
et al. (2006, in preparation) have analyzed the same RGS spec- 
trum of SU Aur using a line-based analysis to derive the EMD 
from the measurement of individual line fluxes. Their results 
are consistent with ours: the EMD peaks at T = 10 MK, with 
an indication of a significant amount of material above ~20 
MK; however, below the peak they find a steeper slope (~ 3) 
than adopted here. 

4.4. Abundances 

The element abundances in the X-ray sources are plotted 
in Fig. |5] The filled circles designate abunda nce ratios rela- 
tive to the solar ph otospheric mix (lAnders & Grevessd ll989: 
iGrevesse & Sau vail 19991 for Fe). The abundances of AB Aur 
show a nearly flat distribution; neither a strong First Ionization 
Potential (FIP) effect (i.e., overabundant low-FIP elements) nor 
a strong inverse FIP effect (i.e., increasing abundances with in- 
creasing FIP) is se en, in contrast to the usual findings in young , 
active stars (e.g.. lArgiroffi et al.ll2004t iTelleschi et alJ 12006). 
Because the emitting material originates, if located in a corona 
or a stellar wind, in the stellar photosphere, a comparison of 
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the X-ray derived abundances with photospheric values will 
be important. Fortunately, a few photospheric abundanc es hav e 
recently been measured for AB Aur dAcke & Waelkensl f2004). 
and the resulti ng normal ized abun dance r atios are shown by the 
open circles 2 . Acke & Waelkens (2004) report a particularly 
low photospheric abun dance of Fe (^ 13% of the solar value 
given by iGrevesse & Sauvallll999l) . The renormalized abun- 
dance distribution is still ambiguous: the coronal abundance of 
Fe is higher than the photospheric value, but a clear FIP depen- 
dent abundance distribution is not visible. 

We also note here that the derived abundances in the SU 
Aur source are atypical for young, active stars and T Tau stars. 
Their X-ra y sources usually show a well-expressed in verse FIP 
effect (e.g.. lTelleschi et alJl2006llArgiroffi etai]l2004l) . In con- 
trast, SU Aur shows O and Ne abundances clearly lower than 
Fe, which is a defining signature of a solar-type FIP effect. And 
second, as suggested earlier, the absolute Fe abundance is quite 
high (0.67-0.81 times the solar photospheric value). Such high 
values are usually reported for relatively inactive stars, while 
magnetic ally active stars reveal muc h stronger depletion of Fe 
llTeUeschi et alll2005UGudelH2004 . 

4.5. Luminosities 

The only previous report on X-ray emission from AB Aur 
that we are aware of is bv lZinnecker & Preibisch (1994). This 
allows us to study possible long-ter m changes in the X-ray 
output. Zinnecker & Preibisch ( 1994) give an X-ray luminos- 
ity of (0.3 ± 0.09) x 10 30 erg s" 1 in the energy range 0.1- 
2.4 keV. In contrast, we find L x « 0.4 x 10 30 erg s" 1 in 
the energy range 0.3-10.0 keV, and L x « (0.5 - 0.6) x 
10 30 erg s~ | in the 0.1-2.4 keV range, t hat is, almost twice 
as much as IZjrme^ker^&^eibischl Jl994l) . We note, however, 
that[zinnecker & Preibisch ( 1994) estimated L x directly from 
count rates, using a conversion factor applicable for a temper- 
ature of 1 keV and an Nh corresponding to Ay = 0.65 mag. 
Modeling the ROSAT count rate for Nh, kT and abundances 
found from our spectral analysis still results in only (0.3 -0.4) x 
10 30 erg s" 1 in the 0.1-2.4 keV range. AB Aur was in a more 
active state during our observation, but we also recall that this 
source is slowly variable on time scales of hours (Fig.^. Our 
luminosity of SU Aur is L x * (7.4 - 7.8) x 10 30 erg s" 1 in the 
0.3-10 keV range. This compares well with values reported by 
ISkinner & Walter I i 19981) . i.e., (8.4 ± 0.09) x 10 30 erg s" 1 in the 
energy range 0.5-10.0 keV. 

4.6. Variability of spectral fit parameters 

In view of the sinusoidal variation of the X-ray count rate in 
Fig. [2 we tested which spectral fit parameters are mainly re- 
sponsible for the modulation. We restricted this study to the 
combined EPIC MOS1+2 data because the RGS signal be- 

2 The photospheric abundances of AB Aur are, with respect to 
the solar photosphere (error ranges in parentheses): N: 1.38; O: 1.62 
(1.38-1.91); Si: 0.78 (0.72-0.83); Fe: 0.13 (0.10 -0.18), referring to the 
solar ab undances of | And ers & G revessd ^989) for N, O, and Si, and 
those of lGrevesse & SauvalHl999l) for Fe. 



comes too weak if the data are split up. We first fitted the entire 
MOS data with two thermal components to obtain a fit very 
similar to the one reported in Table 0] apart from some devia- 
tions in the element abundances (most of which are difficult to 
derive from MOS, in particular O and Ne, given the modest res- 
olution and severe blending). We then split the observation into 
"high state" (first half) and "low state" (second half). Starting 
with the above model fit, we tested whether i) an adjustment 
of Nh (variation due to selective absorption), ii) a renormaliza- 
tion of the emission measures by a common factor (variation of 
Lx), or iii) some changes in all of Nh, kT, and EM are required. 
A variation of Nh could clearly be excluded. This is supported 
by Fig.^that shows that the hard photons vary in concert with 
the soft photons, while they are not significantly affected by the 
weak photoelectric absorption. To quantify this, we have mea- 
sured the amplitude of the sine function relative to the "zero 
level" in each light curve and found them to be the same in 
each energy range, within the errors. No significant changes in 
kT were required, while a simple renormalization of the EMs 
by a common factor produced perfect fits. We conclude that the 
variation of the X-rays are either due to an intrinsic change in 
the luminosity with other plasma parameters remaining equal, 
or due to an energy-independent filtering of photons (e.g., due 
to partial eclipses by the star). 



5. He-like triplets, densities, and the radiation field 

We now discuss the helium-like line triplet of O vn for AB Aur. 
The flux ratio between the forbidden and the intercombination 
lines at 22.1 A and 21.8 A, respectively, is density-sensitive 
roughly in the range of electron densiti es between 10 10 cm" 3 
and 10 12 cirT 3 (Gabriel & Jordan, 1969) for the following rea- 
son: if the electron collision rate is sufficiently high, electrons 
in the upper level of the forbidden transition, ls2s 3 S i, do not 
return to the ground level, Is 2 ^o, instead they are collision- 
ally excited to the upper levels of the intercombination tran- 
sitions, ls2p 3 Pi,2, from where they decay radiatively to the 
ground state. They thus enhance the flux in the intercombina- 
tion line and weaken the flux in the forbidden line. However, 
photons in a UV radiation field may excite the same transition. 
The relevant photon energies correspond to the energy differ- 
ence of the two upper states, and this corresponds to a wave- 
length of 1630 A for the Ovu triplet. The UV radiation field 
is thus important for stars with r e ff of about 10 4 K and more. 
Because T e g of AB Aur has been quoted to be around 10 4 K 
(Table \Tj , we need to co nsider the radiation term. We follow 
iBlumenthal et alJ i 19721) for a rough estimate. The measured 
ratio % = f/i of the forbidden to the intercombination line flux 
can be written as 



1 + 4>/cf) c + rie/Nc i 

where is the limiting flux ratio at low densities CRq « 3.85 
for O vii at the maximum formation temperature), N c is the crit- 
ical density at which H drops to 7?o/2 (N c ~ 3.4 x 10 10 cm" 3 for 
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O vn), and <p c is the critical photoexcitation rate. The influence 
of the radiation field is given by 

cf> _ 3(1 + F)c 3 A(ls2p 3 P 
<p c 8nhv 3 



ls2s 3 S l ) 



(3) 



A(ls2s 3 Si Is 2 l S ) 

where the radiation field energy density is given by Planck's 
equation, 

8tt/iv 3 



W- 



l 



c 3 exp(/iv/£r eff ) - 1 



(4) 



in which we have introduced a geometric dilution factor (e.g. 
iMewe & Schriiverll978l) . 



W 



1 




l/2\ 


1 - 


[-(If 











(5) 



where R is the stellar radius and d is the distance of the source 
from the center of the star (see also lNess et al. 2001 for a more 
detailed presentation). At the surface, W - 0.5. Further, in 
Eq. Q, the A terms are the spontaneous transition probabil- 
ities, (8.12 x 10^ s -1 a nd 1.04 x 10 3 s" 1 , respectively, after 
Blu menfhal et al.lll972l) and F is an expression that we ap- 
proximate by adopting the maximum line formation temper- 
ature for O vn 2 MK ), which yields F = 0.42 for this ion 
(Blu menfhal et all Il972h . We have also assumed that the radi- 
ation field exactly corresponds to the stellar photospheric T e g . 

In a thermal plasma, the flux of the resonance line, r, is 
larger than the flux of the forbidden line, /, and under condi- 
tions relevant to us, namely T > 1 .5 MK, also the sum of / + i 
is smaller than r, G — (f + i)/r < 1 (the "G ratio", see, e.g., 
Pora uet etall koOl). The considerable errors in our measure- 
ments make the extraction of individual lines and the separate 
treatment of "R and Q problematic. To obtain self-consistent K 
and Q while fulfilling other conditions from the spectral fit to 
the entire spectrum, we proceeded as follows. We adopted the 



optimum parameters from the 2-T fit and kept all parameters 
fixed, except for the electron density and (for slight adjustments 
to the total O vn line flux) the emission measure of the cooler 
com ponent. The therm al structure thus sets a requirement on 
Q dPorauet et all 1200 1 ) and also fixes the faint continuum re- 
quired by other lines. We then performed a fit only to the spec- 
tral range around the Ovn triplet, in the wavelength interval 
21.4-22.4 A. 

We first assume negligible influence by the radiation field, 
i.e., we assume that f/i is controlled by the electron den- 
sity. Our fit procedure was making use of the implementation 
of the He-like triplet calculations in XSPEC's vmekal code. 
Varying the electron density changes both the i and / flux 
until a best fit is obtained. We performed this procedure for 
various data binning schemes, using bin widths of 45 mA or 
56 mA. Each time, the best-fit density converged to values 
below the low-density limit for Ovn (n e $ 10 10 cm 3 , cor- 
responding to % — Kq). We then varied n e to find the 68% 
and 90% upper limits. These themselves fluctuated for dif- 
ferent binning schemes; for the average 90% upper limit, we 
find n e ,max,9o ~ (1.3 + 0.4) x 10 11 cm -3 , corresponding to the 
measured f/i = 0.95. For the 68% errors, we find n e ,max,68 ~ 
(4.2 + 1.2) x 10 10 cm" 3 , corresponding to f/i = 2.41. Fig.© 
shows the fit for the low-density limit (solid histogram), and the 
90% upper limit (dotted histogram). The strong / line clearly 
requires a low-density environment; the 90% limit, while for- 
mally acceptable, requires i > /, unlike the data. The densities 
suggested her e, w R < 10 11 cm " 3 , are very typical for stellar 
coronae (e.g.. lNess et all200ll) . 

We now consider the influence of the radiation field. UV 
radiation will lower the // i ratio and thus simulate higher den- 
sities, i.e., the electron densities reported above are overesti- 
mated (Eq. We now assume that the f/i ratio is not sup- 
pressed by high electron densities and ask how far the source 
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Fig. 6. Fit of the Ovu triplet using variable electron density. 
The solid histogram gives the best fit (n e < 10 10 cm -3 ), while 
the dotted histogram is for the 90% upper limit to n e (n e < 
1.3 x 10 11 cm" 3 ). 



must be from the star in order to show the observed f/i ra- 
tio. For r eff = 10050 K (Table we find from Eq. (0-0, 
that the 68% upper limit (f/i = 2.41) is attained at d = 4.6/?, 
whereas the 90% upper limit (f /i = 0.95) requires d > 2.1R . 
For r eff = 9500 K reported bv lvan den Anckeret alJ Jl99«h . 
the 68% upper limit corresponds to d — 3.6R and the 90% up- 
per limit to d — 1 JR. For a given ratio, the origin of most of 
the emission must fulfill these requirements, but some contri- 
butions from closer to the star are not excluded. 

We have assumed here that the radiation temperature is 
identical to 7W which , due t o spectral modifications, may be 
inaccurate. Ness et al. (2001 ) found, for the relevant radiation 
at 1630 A, Trad se veral h undred K below for nearby F and 
G-type stars, while lNess et al J 12002) reported T mA about 300 K 
below r eff for the B8 star Algol (r eff = 13000 K). But even if 
we adopt T mA = 9000 K for AB Aur, we still require d > 2.8R 
and d > 1 3R for the 68% and the 90% limit of the f/i ratio. 

In conclusion, we find that the electron density in the source 
must not exceed a few times 10 10 cm -3 (about 10 11 crrT 3 for 
90% upper limit) and the majority of the source plasma must 
not be closer to the stellar center than (1.3 - 2.1)/?, depending 
on the exact temperature of the radiation field. 



6. Discussion 

In Table H we have listed the basic properties of AB Aur and 
we have compared them with those of the CTTS SU Aur and 
HD 163296, an Herbig star that has been reported to reveal a 
soft spectrum JSwartzet alll2005h . 

The similarity between AB Aur and HD 163296 is striking. 
The stars not only have very similar fundamental properties, 
but their X-ray properties are also very similar. In contrast, SU 
Aur, a relatively massive but markedly later-type classical T 
Tau star, shows X-ray properties very different from the Herbig 



stars, although thes e properties are typical among T Tau stars 
JCiidel et allkoodt . 

We now discuss various models for X-ray production pro- 
posed in the previous literature and compare predictions made 
by these models with ou r spectroscopic results. We l argely fol- 
low the presentation bv lSkinner & Yamauchillll996h . 

6.1. Atmospheric and wind structure 

Numerous optical and ultraviolet spectral observations of AB 
Aur have converged to a wind+chromosphere model in which 
an expanding wind overlies an extended, hot chromospheric 
layer. The latter has a heig ht of about 1.5/?, with a tempera- 
ture pe ak of 1.7 x 10 4 K JBohm & CatalaLll993l) . |Praderie etail 
(1986) and ICatala et alJ (11986) discovered periodic modula- 
tion in the Mg n and O n lines, although the periods disagreed, 
the period of Mgn being 45 hr and the period of On be- 
ing near 32 hr. The latter was identified with the stellar rota- 
tion period, implying formation of O n close to the star, while 
Mgn forms further away from the star (in its wind at sev- 
eral R) and may therefore be modulated by the rotation pe- 
riod of the envelope at that distance. The modulation was 
suggested to be due to a non-axisymmetric wind in which 
fast and slow streams alternate. Magnetic fields would then 
provide a possible explanation to this wind structure. High- 
temperature lmej3_of_Nv anrJO_vi_were_d£tected by FUSE, 
in AB Aur dBouret et all Il997t EToberge et all l200ll) and the 
similar HD 163296 (foeleuil et all 2005). and these were in- 
terpreted as originating fro m shocks fo rmed when the fast and 
slow winds collide (Bou ret et all 1 1997). The same model could 
also produc e X-rays in a layer close to the star (0.05/? above 
the surface, iBouret et all 1 19971) . However, the f/i ratio that 
we measured in the X-ray spectrum suggests that the X-rays 
are formed at distances d > 1.3 R from the center of the star 
(for T rad > 9000 K). Alternatively, Nv and Ovi could also be 
produce d in a wind shock along with X-rays in the model by 
iBabel & Montmerldd 19971) described below. 

Catal a et all i 19991) extended periodicity studies to photo- 
spheric lines and found that the amplitudes of the red emis- 
sion components are modulated with a period of 43 hr, whereas 
the velocity of the blue absorption components is subject to a 
34 hr period. Further, the Hel D3 line shows red and blue com- 
ponents that are both modulated with a period around 45 hr 
(red in amplitude, blue in velocity). These authors interpreted 
the blue Hel D3 components as originating from an equatorial 
wind, whereas the photospheric blue absorption would come 
from high-latitude photospheric regions with radial flows, in- 
dicating a shorter period at high latitudes than near the equa- 
tor. Finally, all redshifted components would be due to polar 
infall. Their 43-45 hr period remains unexplained, but the in- 
fall may be magnetically linked to outflows at lower latitudes. 
Alternatively, the 45 hr Hel (blueshifted) outflow signatures 
may originate from a magnetic disk wind, indicating an anchor 
point of the magnetic fields at 1.6/? (C atala et allll999l) . 

The above interpretation is subject to one caveat, namely 
the assumed high inclination angle (70° in lCatala et al.l fl999). 
Adopting a radius of 2.3 R Q (Table vsin/ = 80 km s _1 
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(Bohm&Catala 1993, Table LL), and the extreme value of 
i = 70° (see also references in ICorder et al.l 12005). we find 
P — 33 hr, in perfect agreement with the 32-34 hr period re- 
ported from O ii and blue photospheric absorption (for i = 90°, 
the maximum rotation period is 35 hr). However, recent ob- 
servations have strongly revised the disk in clination angle and 
now suggest i ~ 21.5° JCorder et alJl2005l and further refer- 
ences therein). This suggests P = 12.9 hr. The 32-34 hr stellar 
rotation period could only be maintained if the disk and stellar 
axis were grossly mis-aligned. 

It is interesting that our X-ray light curve period of 42.5 + 
4 hr perfectly agrees with the period in Mgn and He I, i.e., 
components formed in the chromosphere at the wind base and 
in the wind itself, but is clearly not compatible with the rota- 
tion period based on vsim measurements with ; « 21°, or the 
photospheric blue absorption components. We thus tentatively 
conclude that the X-ray production may be related, in some 
ways, to the wind. 

6.2. Wind shocks 

In hot stars, shocks are driven by instabilities in line-driven 
winds (see, for example. iFeldmeier et aljfl997al for a review). 
The electron temperature of AB Aur's X-ray source is simi- 
lar to those measured in O stars ( Feld meier etallfl997al) . The 
important parameter for the shock instability is the Eddington 
parameter 

t =2J^W m (6) 

llCastoretallll975l) where we set the electron mass scattering 
coefficient cr e = 0.4 cm 2 g -1 , and M is the stellar mass. We thus 
obtain r = 5.6xl0~ 4 «C 1. Under these circumstances, the wind 
cannot be radiation-driven, which makes the instability-shock 
hypothesis unlikely. 

On the other hand, IZinnecker & Preibischl \ 19941) have 
suggested that the wind slams into dense molecular mate- 
rial in the ambient medium. U sing an average wind veloc- 
ity of ICatala & Kunaszl lll987l) . y w = 225 km s" 1 , and the 
wind mass loss rate reported by ISkinner et alJ Jl993l) . M w = 
1.1 x 10~ 8 M Q yr~' from radio observations, we derive a "ki- 
netic luminosity" M w v 2 J2 w 1.8 X 10 32 erg s _1 . That means 
that only about 0.3% of the available wind energy needs to be 
dissipated in shocks to produce the X-ray luminosity. 

However, the observed systematic variability of the X-ray 
source on time scales of several hours makes models based on 
very-large-scale shocks as well as on many distributed shocks 
in the wi nd of a single star u nlikely. This echoes the conclu- 
sions by Skin ner et alJ |2005) who su mmarize observations of 
X-ray variability in O stars. Although Feldmeie r et al J 111997b) 
proposed variability owing to colliding shells in the O-star 
wind, the time scales for such collisions would be shorter (of 
order 500 s). 

6.3. Magnetically confined winds 

Although strong magnetic fields are not expected on Herbig 
stars given their predominantly radiative interior, large-scale 



fossil magnetic fields from the stellar formation process could 
still be trapped in the star. There have been a number of investi- 
gations that studied the consequences of ionized winds trapped 
in l arge-scale stellar magn etic fields. 

lHavnes & Goertzl i ll 9841) presented a model of a global stel- 
lar magnetosphere which is fed by a wind from the stellar 
surface. They considered a magnetic field with a dipole con- 
figuration, where plasma is confined by closed magnetic field 
lines. Close to the stellar surface, where the gravitational forces 
exceed the centrifugal forces, the density is low, while it in- 
crease further out. The temperature also increases, possibly to 
coronal temperatures, at 5-10 stellar radii. The outward trans- 
port of plasma takes place by events where the magnetic lines 
are broken due to excessive density. The energy input in the 
corona from rotation is given by E = 1/2 M W Q?R 2 {L\ - L 2 ), 
where £2 is the angular velocity. This can be rewritten as E « 
1 .7 x \0 11 M W R 2 {L 2 2 - L\) Pf erg s" 1 (where M w is in M Q yr _1 , 
R is in cm and the rotation period of the star P^ is in days). The 
parameters Ly and Li are the distances to the inner and outer 
edges of the X-ray emitting region in units of the stellar radius. 
Both L\ and L2 are a function of the mass loss, the magnetic 
field a nd the rotation per i od, an d are therefore difficult to es- 
timate. lHavnes & Goertzl (11584) estimated these values to be 
L\ = 15 and Li = 20, for a hotter and more massive star. Even 
if we cannot constrain these two values, we expect that (L^-L 2 ) 
will be larger than 1, so that we obtain E > 4 x 10 31 erg s , 
i.e., enough to produc e the observed Ly . 

In another model, llJsov & Melrose! i 19921) considered the 
wind zone (open field lines) outside the corotating magneto- 
sphere (dead-zone). According to their model a current sheet is 
formed in the equatorial plane outside the dead zone that sepa- 
rates regions with opposite directions of the magnetic field. The 
authors estimated the temperature and the energy released in 
the current sheet, assuming bremsstrahlung as a cooling agent. 
Assuming a surface magnetic field strength of 100 G (se e re- 
cent measurements on Herbig stars bv lHubrig et al.l E004). the 
above M w , an average wind velocity of v w = 225 km s _1 
(Catala & Kunasz, 1987; Ski nner et al.lll993l) . and a radius of 
2-4 R , we find for the two parameters £ and 77 in their Eqs. 
(5) and (23), £ » 1 and 77 > 1, r espect ively, and therefore, 
from Eq. (25) in llJsov & MelroseJ (I1992I) . an X-ray luminos- 
ity of 1.1 x 10 30 erg s _1 . Their Eq. (20) predicts an electron 
temperature of T ~ 6 x 10 7 (v lt ,/10 8 cms 4 ) 2 « 3 x 10 6 K. 
These parameters are again in good agreement with our mea- 
surements. 

Finally, Bab el & Montmerlel i 19971) considered the wind 
shock inside the magnetosphere that emerges when the magnet- 
ically guided winds from the northern and the southern hemi- 
spheres collide in the equatorial plane. They predict a shock 
temperature ofT s = 1.13 X 10 5 (v w /10 7 cm s 4 ) 2 * 0.5 - 1 MK 
for v w = 200 - 300 km s _1 (extreme values reported by 
Cata la & Kunaszlll987l) . Magnetically guided winds develop 
shocked equatorial "disks" only if the magnetic fields are suf- 
ficiently strong for confinement. For this, the (equatorial) wind 
confinement parameter, 77 = B 2 R 2 /(Mvw) must be at least unity 
dud-Doula & Owockil200 2: note that B is the surface magnetic 
field). With our stellar parameters (Table [0 and B = 100 G, 
we find 77 « 20, and hence wind shocks can develop. The tem- 
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peratures provided by this model are nevertheless too low with 
respect to the observations. 

In summary, wind-fed magnetospheres may be promising 
to produce the observed X-ray emission although details of 
the magnetic field arrangement and the wind-field interactions 
would need further elaboration for the specific case of AB Aur. 

6.4. Accretion-induced X-rays 

Accretion has recently gained some attention as a possible con- 
tributor to X-ray emission in classical T Tau stars. Accretion 
shocks at the base of magnetic funnel flows may reach high 
temperatures, and high densities (of order 10 13 cm' 3 , llJlrich 



12004 . There has been some evidence for mass inflow in Herbig 
stars from redshifted absorption components in optical lines 



Il976ft . In standard accretion shock models, the shock heats 
up to a few times 10 6 K, with the ensuing X-rays heat- 
ing the underlying photosphere to produc e an UV excess 
JCalvet & Gullbringt Il998i) . iLamzinl ill 999ft concluded that a 
typically small fraction of the X-rays escape from the shock 
that can be seen in soft X-rays. However, it will be im- 
portant to estimate whether the shock is above the photo- 
sphere at all and therefore visible, and this is not normally 
the c ase for T Tau stars with average accretion characteris- 
tics JCalvet & Gullbring i Il998ft . Little variability s hould be 
seen i n Lx or the e l ectron temperature in shocks ( ILamzinl 
1999). Kastner et alJ J2002I) have proposed accretion-induced 
X-ray production for the unusually soft X-ray emission and 
the high densities measured in the spectrum of the classical T 
Tau star TW Hya. In analogy, ISwartz etall ll2005n suggested 
the same scenario for the soft spectrum of the Herbig star 
HD 163296 (Table but high-resolution X-ray spectroscopy 
was not available for this star. 

A further argument in favor of an accretion model was 
brought up bv lStelzer & Schmittl J2004 who argued that TW 
Hya's anomalously high abundances of Ne and N support 
an accretion scenario; depleted metals would condense onto 
grains in the disk, leaving gas enriched in certain other el- 
ements, and this gas would eventually accrete onto the star. 
However, refering to Ne/Fe and N/Fe abundance ratios, there 
are also evolved stars and non-accreting pre-m ain sequence 
stars that reveal high values for these ratios (see lGudelll2 004 
for a r eview). As for the Ne/O abundance ratio. iDrake et alJ 
(2005) indeed found it to be unusually high in TW Hya com- 
pared to other stars. However, such a high Ne/O abundance 
ratio has so far been measured only in TW Hya. In other 
accreting and non-accreting stars, this ratio is found to be 
half as high as in TW Hya (see for example I Argiroffi et all 
l2005t iRobrade & SchmitA 12004 iTelleschi et all l2006h . and 
consistent with va l ues of a large samp l e of late-type stars 
( brake et all 120051: ITelleschi et all Eool . brake et all d2005t) 
suggested that the anomalous Ne/O abundance in TW Hya is 
determined by the higher degree of metal depletion in this older 
star. This ratio can therefore not in general be used as an accre- 
tion signature in young and less evolved TTS. 

What do we know about accretion in Herbig stars? 
Evidence for accretion from a disk is rather indirect. The tem- 
perature of the photosphere that is heated by infalling material 
is very similar to the undisturbed photosphere ( Muzerolle et all 
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irelli et allfl99rilNatta et allkOOOHCatala et allll999h . but 
the interpretation is model-dependent, w ith M possibly e xceed- 
ing 10~ 7 M yr 1 dSorelli et all Il996h . iBlondel et all Jl993l) 
suggested that hydrogen Lya lines in Herbig stars are due to 
infal ling gas , although no signatures were found for AB Aur. 
iGradv etall & 19991) discussed evidence for infalling gas in the 
case of AB Aur. lMuzerolle et al1 lf2004) interpreted B aimer and 
sodium profiles based on magnetospheric accretion models to 
conclude that M * 10~ 8 M Q yr 1 for the Herbig star UX Ori, 
and for a larger sample of stars that excess fluxes in the Balmer 
discontinuity imply M £ 10 -7 M yr -1 . However, the Balmer 
discontin uity excess has been measured to be 0± 0.1 mag in AB 
A ur (lGarrisonHl978t). at b est implying Ma 10~ 8 M yr 1 (af- 
ter Muzerolle et al. 2004). This is supported by measurements 
bvBohm & Catalal dl993^ who find M < 7.5 x 10~ 8 M yr 1 . 
ICatala et all J 1999ft find explicit evidence for near-polar down- 
flows in AB Aur with velocities of about 300 km s , lead- 



Ms y r (see also 



ing to estimates for M of a few times 10 
iBouret & Catalal2000l) . 

The accretion model is specifically favored for AB Aur by 
the measurements of the red component of the He I D3 line by 
ICatala et alJJl999l) . In fact, they found a periodicity in the red- 
shifted line amplitude of 42^1-5 hr, i.e. very close to the period 
that we measure in the X-ray light curve. 

We first need to check whether an accretion rate of order 
M x 10~ 8±1 M yr -1 suffices to explain the observed X-ray 
output if the latter is indeed generated by accretion shocks. 
The accretion luminosity is, ignoring viscous dissipation in the 
disk, L acc = GMM/(2R) * 6 x 10 40 MM/(M e yv^/R erg s" 1 



or L acc 3o ~ 600MM_ 8 /fl where L acc>30 = L acc /(10 M erg s '), 
M = M/Mq, R = R/Rq, and M_ 8 is M in units of 10~ 8 M Q yr -1 
(similar values were reported for HD 163296 given the very 
similar parameters for this star and AB Aur - see Table ^ 
ISwartz et alll2005l) . We conclude that there is sufficient accre- 
tion energy available to produce the X-rays. 

We now estimate what an accretion model would predict for 
X-ray production on AB Aur. The expected shock temperature 
is, from the strong-shock conditions, T = 3v 2 fim p /l6k where 
the upstream flow velocity is approximately equal to the free- 
fall velocity, v ff = (2GM/R) 1 ' 2 (i gnoring centrifugal forces if 
mass is guided along rotating magnetic fields), m p is the proton 
mass, and fi « 0.62 for a fully ionized gas. We thus find 



T a 5.4 x 10 6 ^ [K] 
R 



(7) 



and with the parameters in Tabled T ~ 6 MK, in good agree- 
ment with our measurements. 

To estimate the shock density, we use the strong-shock con- 
dition «2 = 4«j where n\ and «2 are the pre-shock and post- 
shock densities, respectively. We first estimate n \ from the total 
mass accretion rate and the estimated accreting area on the sur- 
face: M a 4nR 2 fvffn e m p where / is the surface filling factor of 
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the accretion flows, or M_g 
find 3 



m 



4x 10 11 M-i 

£3/2^1/2 — 



cm 



1 x lO- n R 3/2 M 1/2 fm. We thus 



(8) 



For the sake of argument, we adopt M-% = 1, which seems 
to be a reasonable value based on the findings summarized 
above, and which is similar to accretion rates of less mas- 
sive T Tau stars. For ty pical filling factors as discussed by 
ICalvet & Gullbrind(ll998l) . i.e., / = 0.1 - 10%, we find n 2 * 
10 12 - 10 14 cm" 3 . 

Although these are densities similar to those measured on 
the T Tau star TW Hya, the O vn triplet we see in AB Aur re- 
quires densities about 100 times smaller. This could only be 
achieved by lowering M to about 10~ 10 M yr~' or by increas- 
ing the accretion area to essentially the entire stellar surface. 
The former possibility is not supported by (at least tentative) 
measurements of limits to M, as summarized above. Accretion 
onto the entire stellar surface is unreasonable given that the star 
accre tes from a disk, and a wind is present (e.g. JPraderie et all 
1986). 

We have not yet considered the radiation field of the A star. 
As shown in Sect. [5] a minimum distance of 1.3 - 2.17? from 
the center of the star is required for the X-ray source to be com- 
patible with the observed Ovn f/i flux ratio. One way out is 
that the shocks are sufficiently shielded from UV radiation. It 
is not clear, then, how X-rays can escape without any absorp- 
tion in addition to the circumstellar photoelectric absorption 
that agrees well with the optical extinction from circumstel- 
lar dust. This absorption is in fact extremely low in AB Aur, 
comp ared to other young stars in the Taurus-Auriga molecular 
cloud dGtidel et alll 2006). 

In summary, then, certain properties of X-ray production 
in accretion shocks may well be explained by simple shock 
models, but there are serious problems with this explanation, 
in particular related to i) the radiation field, ii) the low densi- 
ties measured in the O vn triplet and iii) the lack of any excess 
absorption. There is also little support with regard to selective 
condensation of metals onto grains in the accretion disk. The 
abundance of Fe, thought to be amon g the elem ents that con- 
dense easily onto grains (Stelzer & Schmitt, 2004), is higher in 
the X -ray source than in the photosphere ( Ac ke & Waelkensl 
2004). 



6.5. Coronal X-rays 

Coronae provide the standard explanation for X-rays from stars 
of spectral type F and later, because these stars maintain an 
outer convection zone that drives a dynamo. AB Aur is, how- 
ever, a late-B or early-A type star. 

3 In principle, n 2 could also be derived from the observed L x and 
the shock volume, which derives from / and the shock height, the lat- 
ter be ing dependent on the inflow velocity and the cooling time, see 
ICalvet & Gullbrina |l998). However, such shock models show that 
the X-rays are attenuated by 1-2 or ders of magnitude due to the in- 
falling material, see Lamzin 1 1999). Assuming that all Lx escapes 
would, when scaled with the observed luminosity, lead to an under- 
estimated «2 and an overestimated shock height. 



Models of Herbig stars have indicated that a transient shell 
in which deuterium is burned may develop in Herbig stars, al- 
thoug h this occurs preferent ially in the later-type Herb ig Ae 
stars JPalla & Startled. Il993l) . The calculations of ISiess etaD 
d2000h of pre-main-sequence evolutionary tracks also predict 
the presence of a thin convective layer in young AB stars. Using 
their evolutionary model, we obtain for AB Aur a thin convec- 
tion zone of 0.2% of the stellar radius. 

The situation resembles that of mid-to-late A type stars 
and early-F type stars on the main sequence. Although some 
of these stars are X-ray sources, they are clearly sublumi- 
nous compared to later- type stars (if normalized with L*) de- 
spite their often rapi d rotation, and the X-ray spectra are soft 
(Pan zera et allll999l) . It appears that the thin convection zones 
of these stars are unable to maintain vigorous magnetic dy- 
namos, resulting in soft, solar-like coronae. A similar situation 
may apply to Herbig stars like AB Aur. 

The average temperature of AB Aur is similar to L av 
of moderately active main-sequence s olar analogs (G2-5 V) , 
such as 7T 1 UMa, x 1 Ori, and k 1 Cet jTelleschi et all 12005ft . 
Assuming that AB Aur reveals a similarly structured corona, 
we expect that Lx scales with the surface area. Adopting R = 
2.3R Q for AB Aur, its L x would be * (5 - 6) x 10 29 erg s"\ in 
agreement with the observations. We also note that the densities 
derived from the O vn lines are in good agreemen t with mea- 
surem ents reported for numerous coronal sources jNess et all 
2004). Finally, the modulation of the X-ray light curve ob- 
served in Fig.^could be due to rotation and is thus also con- 
sist ent with the coronal h ypothesis. 

iTout & Pringld (1 19951) have proposed a non-solar dynamo 
that could operate in rapidly rotating A-type stars based on ro- 
tational shear energy. The model predicts, for the time devel- 
opment of the X-ray emission from the associated corona, 



L x (f) = L 



x.o 



t 

1 + - 

to 



-3 



(9) 



where Lx,o is a quantity that depends on stellar mass, radius, 
the amount of differential rotation in the stellar interior, the 
breakup velocity, the coronal heating efficiency, and the effi- 
ciency of magnetic field generation. Follo wing the argumenta- 
tion and the choice of co nstants given in Skin ner et all 12004) 
and lTout & Pringld dl995h we find L x n « 1.6 x 10 31 erg s" 1 . 
For AB Aur, Eq. (3.15) in ITout & Pringle ( 1995J) gives, using 
the same default parameters, to = 6.3 x 10 5 M_^ 2 £ 3 j' 2 _yr =s 
1 .5 x 10 6 yr. Given the star's age of 4 Myr JPeWarf et alll2003t) . 
we expect Lx ~ (3 - 4) x 10 29 erg s . Using our updated 
parameters for AB Aur, we thus confirm the conclusion by 
ISkinner et all J2004I) that this model provides very good agree- 
ment with the observations of AB Aur, although this is not true 
for most other Herbig stars. As the star evolves, the X-ray gen- 
eration would rapidly decay further. 

We found a coronal Fe abundance that is at least equal to the 
photospheric abundance. This is again consistent with a coronal 
model. In more evolved magnetically active stars, the Fe abun- 
dance generally increases toward lower activity and becomes 
comparable w ith the photospheric abundance in inactive stars 
(Gudel l2004 . 
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We note, however, that centrifugal forces exceed gravita- 
tional forces at a distance of 1.68 R from the center of the star 
according to parameters in Table ^ an d assuming an (uncon- 
firmed) rotation period of 12.9 days. For equatorial magnetic 
fields, then, the coronal radius must be less than 1 .68 R, other- 
wise the loops will be unstable JCollier Cameron! Il988l) . This 
condition is only in marginal agreement with the lower limit to 
the source size from the f/i ratio (Sect.|5j unless f rot is larger. 

6.6. X-rays from a companion? 

A majority of stars form in multiple systems, and t his is partic- 
ularly true for Herbig stars (IFeigelson e t all 2003$. A compan- 
ion may therefore also be responsible for the observed X-ray 
emission in AB Aur. Such a companion w ould most likely be a 
T Tau star. Recently. IStelzer et al. (2006a) have studied a sam- 
ple of 17 Herbig Ae/Be stars with Chandra, concluding that at 
least in 7 stars the X-rays could originate from an unresolved 
companion. Only 6 stars are found to be X-ray emitters with no 
visual or spectroscopic detection. Furthermore, the X-ray prop- 
erties in this stellar sample are very similar to X-ray properties 
of CTTS. 

iBehar et alJ |2004) postulated that the X-ray emission de- 
tected from the late-type B star yU Lep in fact originates from an 
unknown pre-main-sequence companion, given the high f/i ra- 
tio measured in its spectrum. For AB Aur, the f/i ratio requires 
the bulk X-ray emission to originate at r > 1.38/?. This would 
be fulfilled if the X-ray emission originated from a companion. 

Stringent constraints have be en discussed in the literature, 
interpreted and summarized by Piet u et all J2005I) (see their 
Sect. 5.1, and references therein): Any co-eval companion 
within 120-1500 AU (0.86-12.5") must have a mass < O.O2M . 
The mass upper limit is O.25M down to a separation of 0.4". 
In the range of 0.07-10", the upper limit to a companion 
mass is, from spec kle interferometry, 0.05-0. 3M Q . Recently, 
iBaines et al.l J2006I) reported evidence for binarity of several 
Herbig stars, including AB Aur. However, they estimated a sep- 
aration of 0.5-3.0" for AB Aur. 

The point-spread function of XMM-Newton does not allow 
us to distinguish between our target source and potential com- 
panions within a few arcsec. We have therefore analyzed an 
exposure of the region around SU Aur + AB Aur obtained 
from the archive of the Chandra X-Ray Observatory, revealing 
much better positional information (obs ID = 3755). Standard 
data reduction methods and up-to-date aspect corrections were 
applied. The observation used the ACIS-S detector, with the 
high-energy grating inserted. It was centered on SU Aur, with 
AB Aur being located close to the chip edge. Using the wavde- 
tect (wavelet detection) routine in the Chandra CIAO software, 
we measured the centroid positions of the images of both stars, 
to find 



AB Aur 



SU Aur 



RA(J2000.0) 
5(J2000.0) 

RA(J2000.0) 
5(J2000.0) 



= 4 h 55 m 45.846 s 

= 30° 33' 04.138" 

= 4 h 55 m 59.389 s 

= 30° 34' 01.297" 



+0.17" 
+0.17" 

±0.11" 
+0.11" 



AB Aur 



SU Aur 



RA(J2000.0) 
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4 h 55 m 45.826 s 
30° 33' 04.37" 



RA(J2000.0) = 4 h 55 m 59.381 s 
5(J2000.0) = 30° 34' 01.56" 



+0.07" 
+0.07" 

+0.06" 
±0.06" 



The position offsets of the X-ray sources are thus 0.258" and 
0.103" in RA for AB Aur and SU Aur, respectively, and 
-0.232" and -0.262" in declination. The similar offsets for 
both stars suggest a systematic offset of the pointing of order 
0.28", well within the errors of the Chandra attitude solution. 4 
Correcting the AB Aur position by the offset of the better deter- 
mined SU Aur position, the deviation from 2MASS of AB Aur 
is only 0.16" in RA and 0.03" in declination, which is within 
the errors of the measurement. We conclude that we have iden- 
tified AB Aur in X-rays well within 0.5" (at the 3cr level) of the 
2MASS position. It is therefore improbable tha t the source o f 
the X-rays is the companion detected bv lBaines et alJ J2006), 
which is thought to have a separation with AB Aur of > 0.5". 

Further arguments favor intrinsic X-ray emission from the 
Herbig star. In particular, the observed low average temperature 
is rather uncommon to lower-mass T Tau stars; the latter rather 
show hot components with characteristic temperatures up to 
20-30 MK as, for example, SU Aur b ut also other CTTS in the 
XEST survey jTelleschietalll2006h . Very-low mass stars and 
brown dwarfs do reveal softer spectra with dominant tempera- 
tures below 10 MK but their total X-ray lu minosities are much 
below Lx measured here for AB Aur (see iGrosso et alj f2006 
for the brown dwarf sample from the XEST survey). 

We also note that the Nh value d etermined in the XEST 
survey (Nh ~ 5 x 10 20 cm' 1 . iGiidel et alJl2006l) is in perfect 
agreement with the v isual extinction measured for AB Aur by 
iRoberee et alJ J200ll) (0.25 mag), if we apply a st andard inter- 
stella r conversion law, Nh ~ 2 x 10 21 Ay mag JVuong et alJ 
120031 and references therein). 

The most substantial argument against the companion hy- 
pothesis is the close coincidence between the X-ray period and 
the period observed in the lines of the wind of the Herbig star. 

Finally, AB Aur and HD 163296 are very similar in both 
their intrinsic properties and their X-ray properties. If the X- 
rays would indeed originate from nearby T Tauri companions, 
then the companion of AB Aur would happen to be very similar 
to the companion of HD 163296. 

Taken together, these arguments suggest that the X-rays are 
not originating from a companion but from AB Aur itself. This 
is different from flaring, hard sources among Herbig stars for 
which T Tau companions have recently been identified (see 
Sect.^. We therefore suggest that the unusually soft emission 
is indeed a distinguishing property of genuine Herbig star X- 
ray emission. 

6.7. Disk related models 

The X-ray luminosities of A and B st ars decrease as they ap- 
proach and reach the main sequence ( Stel zer et al.L l2006al) : at 
the same time the outflow activity is believed to cease and the 



The nearest 2MASS objects are located at 



4 The 90% source location error circle in Chandra has a radius of 
about 0.5", see Chandra Proposers' Observatory Guide v.8 
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dense surrounding gas dissipates. That suggests that the X-rays 
of AB Aur could be related to the presence of the circumstellar 
disk. Two different models that link the X-ray emission with 
the presence of the disk have been discussed in the literature: 
a disk corona and reconnection of magnetic fields that link the 
star to the disk. We describe these models, although they do not 
currently make predictions that we can test against our data. 

The presence of a di sk corona was discussed by 
Izinnecker & Preibischl i 19941) . but very little is known about 
its generation. The ionization due to the decay of radioactive 
nuclides could increase the conductivity of the disk to a level 
high enough to generate a magnetic field. The differential rota- 
tion in the disk could then generate a disk corona. 

X-ray emission generated by reconnection of magnetic 
fields linking the star to the disk has been discussed as a model 
for the X-ray generatio n mechanism of low-mass protostars 
(Mo ntmerle et all l2000t) . If the rotation period of the star is 
not the same as the rotation period of the disk, a magnetic loop 
connecting the disk with the star will twist and inflate until it 
comes into contact with itself and reconnects. A similar model 
could apply to Herbig stars, although we would expect higher, 
flare-like temperatures due to the reconnection process. 



7. Conclusions 

We have presented the first high-resolution X-ray spectrum 
of an Herbig Ae/Be star, namely AB Aur. The use of high- 
resolution spectroscopy has allowed us to obtain important 
spectral information that cannot be addressed with EPIC spec- 
tra alone. The O vn triplet constrains the electron density and 
is therefore important for the discussion of the different mod- 
els. Further, we have been able to reliably determine the abun- 
dances of the high-FIP elements O and N. Finally, the Ovm, 
O vn and N vn lines permitted us to constrain the cool plasma. 

We found the X-ray spectrum of AB Aur to be rather soft, 
with spectral-fit results that are consistent with a mean coronal 
temperature of about 5 MK, i.e., much less than the usual tem- 
peratures of coronae of low-mass pre-main-sequence stars that 
usually exceed 10 MK. We found an X-ray luminosity of about 
4 x 10 29 erg/s in the 0.3-10 keV range. We derived the abun- 
dances and found them not to follow a First Ionization Potential 
(FIP) distribution, nor an inverse FIP distribution. We normal- 
ized the coronal abundances to the new p hotospheric abun- 
dances of AB Aur found by lAcke & Waelkensl d2004h . who 
measured a very low photospheric Fe abundance. The Fe coro- 
nal abundance then is at least as large as the photospheric value. 

The density-sensitive O vn triplet has been studied in detail. 
Although its S/N is moderate, we found that the line flux ratios 
indicate densities n e 5= 10 11 cm -3 , with the best-fit value being 
at the low-density limit (n e » 10 10 cm"" 3 ), similar to what is 
commonly found in stellar coronae dNess et al l l2004 . 

We have discussed several X-ray generation mechanisms, 
and provided supporting evidence or pointed at problematic 
features for each of them. First, the probability that the X-ray 
emission is originating from a companion TTS is small. The X- 
ray source is identified within 0.5" of the 2MASS source cor- 
responding to AB Aur in the Chandra image, so that the com- 



panion would have a mass < 0.3M o according to the constraints 
given by Pietu et al. (2005). Such a low-mass star would rarely 
produce an X-ray luminosity as high as observed for AB Aur. 
Furthermore, the close coincidence of the period that we have 
measured in X-rays with the period measured in lines formed 
in the wind of AB Aur itself, makes the companion hypothesis 
unprobable. 

Accretion-induced X-ray emission has been widely dis- 
cussed in the recent literature. With the observed electron den- 
sity, this would be possible only if M m 10~ 10 M yr~', 
i.e. lower than the value suggested in the literature (M « 
10~ 8±1 M Q yr _1 ), and a filling factor ^ 10%. The major problem 
with this scenario is that the UV radiation field (of the stellar 
photosphere and the shock itself) would suppress the forbidden 
line of the O vn triplet. 

The hypothesis that X-ray emission is generated by shocks 
in a line-driven wind, similarly to the mechanism that is be- 
lieved to produce X-rays in O stars, is ruled out by the observed 
variability and the inability of the radiation field to drive the 
wind. 

Magnetic fields have recen t ly be e n detected on sev- 
eral Herbig stars dPonati et all 1 19971 iHubrig et all 120041 
IWade et alf E)05). Two further possible mechanisms are fun- 
damentally dependent on the existence of a magnetic field: 
coronal emission and magnetically confined winds. 

Coronal emission of the type seen in the Sun requires the 
presen ce of a dynamo. The recent calculations bv ISiess et all 
(2000) predict a thin convective layer for Herbig stars, which in 
the case of AB Aur amounts to 0.2% of the radius. The question 
then is whether this thin convection layer would be sufficient to 
generate the dynamo that results in a corona with the observed 
X-ray properties. The corona should be quite extended in order 
to allow the f/i flux ratio in the Ovu triplet to be larger than 
unity. This is plausible because the surface gravity of AB Aur 
is about half that of the Sun. 

Magnetically confined winds are a promising alterna- 
tive the more so t hat winds have explicitly been measured 
(Cat ala et all ll999). There is interesting evidence supporting a 
model of this kind. The period of the modulation that we mea- 
sured in the X-ray light curve (42.2 hr, Sect. 14. 1> is very close 
to a modulati on period of Mgn lines thought to form in the 
wind (Praderi e et all fl9*8 6). Whatever the production mecha- 
nism of the X-rays, it is therefore suggestive to assume that 
they are closely related to the wind or are produced by it. The 
advantage of such a model is that X-rays are produced at some 
distance from the stellar surface, which alleviates the problem 
with the suppression of the O vn / line flux. 

In both these last two models, we have encountered a 
problem when we have adopted a stellar inclination angle 
of i » 21.5° dCorder et all 120051) and vsiru « 80 km 
s -1 (Bohm & Catala, 1993). The rotation period would then 
be 12.9 hr. This value is significantly smaller than both the 
Mg n and X - ray pe riod (« 42 hr) and the value suggested by 
ICatala et all |l999) for the rotation period (32-34 hr). Further, 
in the hypothesis of a stellar corona, the latter should in this 
case not exceed a height of 0.68 R (the location of the co- 
rotation radius where centrifugal forces Q.d 2 equal gravita- 
tional forces GM/d 2 ) above the surface at the equator, be- 
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cause otherwise c entrifugal forces would make loops unstable 
(Col lier Cameronl 11988). Further, with the adopted small in- 
clination angle, no significant rotational modulation should oc- 
cur, and second, the f/i ratio is at risk of being significantly 
suppressed. 

The problems with a wind in a magnetosphere or a corona 
would be alleviated if the stellar rotation axis were, for some 
reason, inclined against the axis of the circumstellar disk, or 
the latter were largely warped between the inner r egions and 
the well-observed outer regions (Corder et all 120051) . The 33 hr 
period suggested by ICatala et al.l (119991) to be identified with 
the rotation period would then require a stellar inclination of 
about 70 deg, making partial eclipses of X-ray emitting mate- 
rial, located for example close to the equatorial plane of the star, 
easily possible. At the same time, azimuthally varying wind ve- 
locities would prod uce the line shift periodicity as reported in 
the optical and UV (Praderie et all Il986t ICatala et all fl986). 
We have, however, no explanation as to why the rotation axis 
should be inclined against the disk. 

We further note that HD 163296 displays very similar prop- 
erties to AB Aur (Table [0, with a soft X-ray spectrum simi- 
lar to the one that we have found for AB Aur. However, the 
X-ray properties found for AB Aur and HD 163296 are not 
common to all Herbig stars. Some of these stars in fact dis- 
play harder spectra, with temperatures rea ching several tens 
of MK and with larger X-ray luminosities dHamaguchi et all 
l2005USkinner et all l20o4 IStelzer et alll2006ah . In fact, given 
the peculiar properties of the X-rays of AB A ur and the sim- 
ilar HD 163296 (Table □ Swartz et al.ll2005l) . namely a very 
soft X-ray spectrum with a moderate Lx, we suggest that these 
properties are defining properties of X-ray emission intrinsic to 
Herbig stars, while hard, flaring emission may be due to unde- 
tected companions (see Sect.^ an d Sect. 16.61 . 
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